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ABSTRACT 



o 
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o 



We analyse visual observations of the pulsations of the red giant variable L2 Pup. The 
data cover 77 years between 1927 and 2005, thus providing an extensive empirical base for 
characterizing properties of the oscillations. The power spectrum of the light curve shows a 
single mode resolved into multiple peaks under a narrow envelope. We argue that this results 
from stochastic excitation, as seen in solar oscillations. The random fluctuations in phase also 
support this idea. A comparison with X Cam, a true Mira star with the same pulsation pe- 
riod, and W Cyg, a true semiregular star, illustrates the basic differences in phase behaviours. 
The Mira shows very stable phase, consistent with excitation by the K-mechanism, whereas 
W Cyg shows large phase fluctuations that imply stochastic excitation. We find L2 Pup to 
be intermediate, implying that both mechanisms play a role in its pulsation. Finally, we also 
checked the presence of low-dimensional chaos and could safely exclude it. 
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1 INTRODUCTION 

L2 Puppis (HR 2748; HIP 34922) is a bright nearby red giant with a 
pulsation period of about 140 d. In Paper I (Bedding et al. 2002), we 
showed that the star has undergone a remarkable change in mean 
visual magnitude over the past century, and is currently undergoing 
a dramatic decline. We argued that the most likely cause is the re- 
cent formation of dust in an extended atmosphere. In this paper, we 
discuss the pulsation behaviour of the star and suggest that convec- 
tion makes a significant contribution to the driving mechanism. 

The mechanism by which pulsations are excited in long- 
period variables is a subject of current interest. Mira variables have 
large amplitudes and are very regular, reflecting the nature of the 
driving process, which is self-excitation via opacity variations (i.e. 
the K-mechanism). Semiregular variables (SRs), by contrast, have 
smaller amplitudes and less regular light curves. The nature of these 
irregularities is far from being understood. Multiple periodicity has 
been found in many cases and is generally interpreted in terms of 
multimode pulsations (Kiss et al. 1999, Percy et al. 2001). This in- 
terpretation was also supported by the multiple period-luminosity 
relations of SRs in the Magellanic Clouds (Wood et al. 1999, Wood 
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2000) and by observed mode changes in several semiregular stars 
(Cadmus et al. 1991, Percy & Desjardnis 1996, Bedding et al. 1998, 
Kiss et al. 2000). However, even the clearest examples of multiply 
periodic SRs exhibit seemingly irregular additional changes in the 
light curves (Lebzelter & Kiss 2001, Kerschbaum et al. 2001), so 
that multimode pulsation with stationary frequency and amplitude 
content does not fully explain the light variations in SRs. Other ex- 
planations suggested in the literature include chaotic phenomena 
(Icke et al. 1992, Buchler et al. 2004), coupling of rotation and 
pulsation (Barnbaum et al. 1995, Soszynski et al. 2004) and dust- 
shell dynamics (Hiifner et al. 1995, 2003). Studies of radial velocity 
variability in SRs (Lebzelter et al. 2000, Lebzelter & Hinkle 2002) 
suggested that light variations are dominated by stellar pulsations, 
which implies that the observed irregularities must be largely re- 
lated to stochastic behaviour of the oscillations. 

In semiregulars, it seems plausible that there is a substantial 
contribution to the excitation and damping from convection. In- 
deed, Christensen-Dalsgaard et al. (2001) have suggested that the 
amplitude variability seen in these stars is consistent with the pulsa- 
tions being solar-like, i.e., stochastically excited by convection (see 
Bouchy & Carrier 2003 and Bedding & Kjeldsen 2003 for recent 
reviews of solar-like oscillations). Note that we are following con- 
vention by referring to oscillations as 'solar-like' if, like those in 
the Sun, they are stochastically excited and damped. This is in con- 
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Lj Pup (1927 -2005) 




L Pup (JD=2425000 + n ■ 2375 ) 



Figure 1. The full light curve of L2 Pup, showing the gradual mean bright- 
ness variations as 10-day means (top) and the coiTesponding amplitude 
spectrum (bottom). The small insert shows the spectral window (in am- 
plitude), with sidelobes at it 1 cycle per year (0.0027 c/d) reflecting annual 
gaps in the data. 



trast to oscillations in "classical" pulsating stars, such as Miras and 
Cepheids, in which the excitation occurs via a feedback in opacity 
variations (the k mechanism). We stress that in other ways, the os- 
cillations in red giants are definitely not solar-like. For example, the 
Sun oscillates in many modes, both radial and non-radial, whereas 
red giants oscillate in only a handful of modes that are thought 
to be exclusively radial. In K giants, the subject of Mira-like ver- 
sus solar-like excitation of the oscillations has been discussed by 
Dziembowski et al. (2001). They attempted theoretical modelling 
of the oscillations in a UMa and other red giants and suggested that 
there might even be a mixture of Mira-like and solar-like modes 
in certain cases. The question of solar-like oscillations in SRs was 
discussed by Bedding (2003), who interpreted the Lorentzian shape 
of the power spectra of several bright SRs as evidence for stochas- 
tically excited oscillations. Here we further develop this idea by 
analysing the decades-long visual light curve of L2 Pup. 



2 DATA ANALYSIS 

The updated light curve of L2 Pup is plotted in the top panel of 
Fig. Q The observations since JD 2451000 were made by A.J. 
Jones, P. Williams and L.L. Kiss. The dataset consists of 1981 
points between JD 2425249 and 2453487, each one being a 10- 
day average. This is roughly 1000 days longer than in fig. 1 of 
Paper I and it clearly shows that the latest dimming event is still 
in progress, though it may have slowed down a little bit. The fre- 
quency spectrum (bottom panel in Fig. is dominated by a set 
of low-frequency peaks, caused by the slow variations of the mean 
brightness, and the main group of closely separated peaks at the 
pulsation frequency ~0.007 c/d. There is also a slight indication 
of an increased amplitude density around ~0.015 c/d, which may 
be attributed to either a harmonic of the main pulsation (indicating 
non-sinusoidal variations) or a second (overtone) pulsation mode. 

Details of the light curve can be seen in Fig.|2| Here, we sub- 
tracted the mean trend of the light curve by fitting a low-order poly- 
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Figure 2. Oscillations in L2 Pup (after removal of the long-term trend). 
Each segment shows a 2375-d long subset, which covers roughly 17 cycles. 
The mean photometric uncertainty is about the same as the symbol size. 



nomial to the data. Individual cycles have peak-to-peak amplitudes 
ranging from practically zero up to over 2 magnitudes, with no cy- 
cles having the same shape. This illustrates very well the semiregu- 
lar nature of L2 Pup. There seems to be no correlation between the 
actual mean brightness and amplitude or cycle shape, which sup- 
ports the idea that gradual brightness changes are not coupled to 
the pulsations. 

The data were submitted to the following analysis: 

( i) the mode lifetime was determined from a power spectrum fit; 

(ii) phase variations were compared to a true Mira star with the 
same period; 

( Hi) the presence of low-dimensional chaos was tested with phase 
space reconstruction. 



3 DISCUSSION 

3.1 Power spectrum and mode lifetime 

The power spectrum of the light curve is shown in Fig.|3| We show 
a close-up of the region corresponding to the 140-day pulsation 
period, and see that the power is split into a series of peaks under a 
narrow envelope. Note that the long term trends in the light curve, 
which we discussed in detail in Paper I, are irrelevant here because 
they occur at much lower temporal frequencies. We have verified 
that removing these long-term variations by subtracting a low-order 
polynomial has neglible effect on the power spectrum in Fig.|3| 

The power distribution in Fig.|3|is typical of a stochastically 
excited oscillator and is strikingly similar to close-up views of 
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Figure 3. Close-up view of the power spectrum around the frequency of 
pulsations. 



individual peaks in the power spectrum of solar oscillations. As- 
suming a stochastically excited damped oscillator, we have fitted a 
Lorentzian profile assuming statistics with two degrees of free- 
dom. This is based on a maximum likelihood fit assuming an expo- 
nential distribution of the noise (Anderson et al. 1990, Toutain & 
Frohlich 1992). The fit gives a period of 138.3 d with a half-width 
at half power of F = 9.0 x 10"^ d^^, which implies a mode life- 
time of r = (27rr)~^ = 4.8 yr (i.e. approximately 12.5 pulsation 
cycles). 

This represents one of the first measurements of mode lifetime 
(or damping time) in a star with solar-like oscillations. Apart from 
the Sun itself, in which the lifetime is a 2-4 days (e.g. Chaplin et 
al. 1997), there are recent measurement of oscillation lifetimes in 
four stars. Bedding et al. (2004) and Kjeldsen et al. (2005) reported 
mode lifetimes of 2-3 d in both a Cen A and B. For the G giant 
5 Hya Stello et al. (2004) suggested a lifetime of ~2 days (about 
20 mean cycle lengths), which is much shorter than was calculated 
theoretically (Houdek & Gough 2002). Finally, in the K giant Arc- 
turus, which is the closest example to L2 Pup, photometric data 
taken by the star tracker on the WIRE satellite implied strongly 
damped oscillations with a lifetime comparable to the mean period 
of 2.8 days (Retter et al. 2003). Since the damping rate depends 
on the stellar structure and convection properties in a complex way 
with a number of weakly constrained theoretical parameters (e.g. 
Balmforth 1992), and moreover, there are no theoretical calcula- 
tions directly applicable to L2 Pup, it is currently impossible to 
qualify the agreement with theoretical expectations. We will ad- 
dress the dependence of mode lifetime on main physical parameters 
in a subsequent study of bright semiregular variables. 



3.2 Oscillation amplitude and phase 

Could the oscillations in L2 Pup be driven entirely by convection, 
without any input from the «:-mechanism? This has been suggested 
for semiregular variables by Christensen-Dalsgaard et al. (2001), 
who described such oscillations as solar-like. The expected ampli- 
tude is highly uncertain, but we may extrapolate from solar-type 
stars using empirical scaling laws. Visual photometric amplitudes 
of red giant stars are highly dependent on the mean temperature, 
whereas velocity measurements, when available, give a more direct 
measure of the physical pulsation amplitude. Kjeldsen & Bedding 
(1995) proposed that the velocity amplitude of solar-like oscilla- 




Figure 4. Phase variations in L2 Pup, W Cyg and X Cam. C(t) and S{t) 
are the components of the Gabor transform, which give the instantaneous 
phase as ip{t) = arctan(5(t)/C(t)). 



tions should scale as L/M. More recently, Kjeldsen & Bedding 
(2001) suggested a revised scaling relation (to account for low ob- 
served amplitudes in F stars), in which velocity amplitude scales as 
l/g (which is proportional to I//(Mr*g )). 

Adopting L = 1500 Lq, T^a = 3400 K and M = 1 Mq 
(Jura et al. 2002), the original and revised scaling relations pre- 
dict velocity semi-amplitudes for L2 Pup of about 400 m/s and 
3 km s^^, respectively. From the observational side, there are two 
published radial velocity measurements for L2 Pup: Cummings et 
al. (1999) measured a semi-amplitude of about 2.5 km s^^, while 
Lebzelter et al. (2005) reported six data points in a full range of 
12 km s~^ over almost three pulsation cycles, implying a semi- 
amplitude of 6 km s~^. Given the large uncertainties of the extrap- 
olated relations, the observed amplitudes do not contradict predic- 
tions. It therefore seems possible that the oscillations in L2 Pup are 
excited entirely by convection, with no Mira-like contribution. 

We find further evidence for non-Mira behaviour from a com- 
parison of phase changes in L2 Pup with those of a true Mira with 
a very similar pulsation period. A well-observed northern circum- 
polar star is X Cam, for which the General Catalogue of Vari- 
able Stars (Kholopov et al. 1985-1988) fists my = 7™4 - 14™2, 
P=143.56 d. We downloaded its visual data collected by the Asso- 
ciation Frangaise des Observateurs d'Etoiles Variables (AFOEV^), 
extending back to the early 20th century, and binned them into 10- 
day averages. Also, as an example of a "real" semiregular, we took 
visual data of W Cyg, an SRb type star with two dominant peri- 
odicities, 131 d and 235 d (Howarth 1991). Then we calculated the 
following quantities for each star; 



C{t) = ^w{t,ti){m{ti) - {m(U)))cos{2-Kfti) 



(1) 
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Figure 5. Phase vaiiations in L2 Pup, W Cyg and X Cam as function of 
time. 
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which are closely related to the Fourier transform of the light curve 
{m{ti)} (i = l...n). The only difference is the presence of the 
Gaussian weight-function w{t,ti), which was used as a moving 
window over the light curve (in other words, C{t) and S{t) are the 
real and the imaginary components of the Gabor transform; for a 
recent review of related transforms see Buchler & KoUath 2001). 
In our case, w{t, ti) was centered at t and ran between t\ and tn 
with a time step of 50 days, and the full-width at half maximum 
was fixed as 700 days (about 5 cycles). This way we could mea- 
sure the local phase tp{t) = arctan(5'(f)/C(f)) at fixed frequency 
f — 1 /P. For X Cam we determined P=143.69 d, very close to the 
GCVS value, while for L2 Pup we used the mean period given by 
the center of the fitted Lorentzian (138.3 d). In the case of W Cyg, 
we took the shorter period (131 d), which is the dominant one. As 
usual, {m{ti)) denotes the mean magnitude. 

First we show the phase variations by plotting S{t) versus 
C{t) in Fig.|4| The fact that L2 Pup and W Cyg are closer to the 
origin than X Cam is a simple consequence of their smaller light 
curve amplitudes (proportional to y^C{ty + S{ty). The interest- 
ing thing is the phase range covered by the stars. While X Cam was 
meandering over a fairly narrow range of about 20-30°, drawing 
a loosely defined arc in the phase plane, L2 Pup went almost all 
around the full circle. We interpret this difference as evidence that 
in L2 Pup, the changes in both amplitude and phase are close to 
random, almost as randomly as in W Cyg, which is centered on the 
origin. 

The intermediate nature of L2 Pup is also apparent from the 
phase variations as functions of time (Fig.|5j. X Cam's phase is 
very stable in time, with no sudden changes; L2 Pup has good 
phase stability over several thousand days but then some jumps oc- 
cur. In contrast, W Cyg shows random phase variations on short 
time scales (see Howarth 1991 for further discussion of the phase 
changes in this star). The phase variations of X Cam are fully 
consistent with the K-mechanism, which by definition includes a 
phase coherent positive feedback from periodic opacity changes. 
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Figure 6. A typical light curve segment with the smoothed and interpolated 
signal. 



The phase of W Cyg varies continuously, supporting the view that 
the oscillations are stochastically excited, presumably by convec- 
tion. Finally, L2 Pupis an intermediate case. The phase fluctuations 
imply stochastic behaviour, but there may also be some driving 
from the K-mechanism. 



3.3 A test for low-dimensional chaos 

Could the seemingly complex pulsational behaviour of L2 Pup 
be caused by a simple low-dimensional chaotic system? Recently, 
Buchler et al. (2004) concluded that the irregular pulsations in 
some semiregulars are the result of the non-linear interaction of two 
strongly nonadiabatic pulsation modes, although they admitted that 
stochastic procedures may also affect the light curve shape to some 
extent. The first clear detection of chaos in a Mira star was pre- 
sented by Kiss & Szatmary (2002), who demonstrated the chaotic 
origin of the amplitude modulation of R Cygni. Here we employ 
the same non-linear tools of light curve analysis in order to yield 
insights into the dynamics of the pulsations in L2 Pup. Instead of 
using the whole dataset, we restricted the nonlinear analysis to the 
better-sampled two-thirds of the data, from JD 2435000 onwards. 

We made the following pre-processing steps: 1. long-term 
trend removal by subtracting a polynomial that was fitted to the 
binned light curve; 2. smoothing and interpolating the resulting 
light curve. The first step was made to remove the long-term mean 
brightness variations that are caused by the increased dust extinc- 
tion events and are not connected with pulsation (Paper I). We ex- 
perimented with different orders of polynomials and accepted 10 
as the best one which followed quite well the sudden changes of 
the mean level. After subtracting the polynomial, the residual light 
curve clearly showed the complex amplitude changes. Further in- 
terpolating consisted of fitting Akima splines to the data with the 
constant sampling rate of 10 days. The spline fit was smoothed with 
a Gaussian weight-function with an FWHM of 20 days. We show a 
part of the final result in Fig.|6| 

The pre-processed light curve was submitted to our nonlinear 
analysis. The applied numerical tools and their basics have been 
described by Hegger et al. (1999). Briefly, we used time-delay em- 
bedding to reconstruct the phase space of the system. First, we es- 
timated the embedding parameters, the optimal time-delay embed- 
ding dimension. We then calculated various projections to visual- 
ize the trajectories in the phase space. The existence of an attractor 
would produce regular structures, and finding these was the aim of 
the analysis. 

As for R Cygni, the optimal delay was found to be between 
10-30% of the formal period. Since the reconstructed images did 
not show strong dependence on the time delay, we adopted 50 days. 
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Figure 7. Broomhead-King projections. See text for further details. 



corresponding to 5 steps in the dataset and about 30% of the mean 
period. We tried different embedding dimensions (de) and the most 
informative images were found for de=4. To visuahze the embed- 
ded state vectors, we have generated the Broomhead-King projec- 
tions (Broomhead & King 1986), which project onto the eigenvec- 
tors of the correlation matrix. 

The most informative projections are shown in the top row of 
Fig. For comparison, we also show the same projections for a 
shorter dataset without long-term trend removal but with the same 
smoothing and interpolation procedure. This subset was chosen be- 
tween JD 2440000 and JD 2450000, when the mean brightness of 
L2 Pup was fairly constant. Both datasets yield the same result: 
there are no regular structures in the reconstructed phase space of 
L2 Pup and the presence of an attractor can be safely excluded. The 
system is temporarily close to stable states from time to time (as 
suggested by the circular structures in the first column of Fig. 0, 
but it then changes erraticly and the phase space projections be- 
come fuzzy clouds. We have checked this conclusion up to dc — 10 
and the results did not change, producing only similar blurred struc- 
tures in higher dimensions. We conclude that the complex light 
variations cannot be explained by low-dimensional chaos. 

We also note that Buchler et al. (2004) made a somewhat ex- 
trapolative argument against stochasticity in semiregular variables. 
While criticizing Konig et al. (1999) for interpreting R Scuti's 
light curve in terms of stochastic damped oscillators, Buchler et al. 
(2004) rejected the idea as being meaningless on physical grounds, 
because no mechanism is proposed that could excite damped modes 
to such large amplitudes, like the factor of 40 in the light curve 
of R Set (i.e. that star's visual amplitude can reach 4 mags). Fur- 
thermore, they also considered the energetics of a 0.7 M0 pulsa- 
tion model with L=1000 Lq and Tcff=5300 K (implying 37 Rq for 
model radius). Their main argument was the fact that the model's 
average pulsational kinetic energies in alternating large and small 
amplitudes exceeded by a factor of two the average turbulent ener- 
gies, so that even if all available turbulent energy were converted 
into pulsational kinetic energy, there is not enough energy. How- 
ever, while this consideration might be true for R Set, which is an 
RV Tauri type pulsating yellow supergiant, both the large amplitude 
in R Set and the quoted model properties are highly irrelevant when 
considering pulsations in red giant stars. Firstly, the visual ampli- 



tude of semiregular and Mira stars are grossly enhanced by the ex- 
treme temperature sensitivity of molecular species in their much 
cooler atmospheres with typical temperatures around 3,000-3,300 
K. In addition, the visual range is in the Wien-part of the spectrum 
for these stars, which also increases the amplitude's temperature 
sensitivity. For that reasons, infrared observations (e.g. in K band) 
draw a more reliable picture of their luminosity fluctuations. Typ- 
ical near-infrared amplitudes in semiregular stars are between 0.1- 
0.3 mag (Whitelock et al. 2000, Smith 2003), which means there is 
no need for any exotic mechanism to excite damped modes to large 
amplitudes. For L2 Pup, the DIRBE instrument on the COBE satel- 
lite measured less than 10% fluctuations at 2.2/i (Jura et al. 2002, 
Schmidt 2003), which is typical among Mira-like semiregulars (i.e. 
SRa type SRs). Secondly, we feel that the applicability of a 5300 K 
model to 3000-3300 K stars must be limited and speculative. Red 
giant stars can have luminosities and radii larger by up to a factor of 
ten, in which parameter region the properties of convection are not 
necessarily the same as in hotter and smaller stars. Consequently, 
while Buchler et al. (2004) presented convincing evidence for low- 
dimensional chaos in three (or maybe four) semiregular stars, their 
arguments against stochastic excitation do not lie on firm physical 
basis. We therefore believe that stochastic (solar-like) oscillations 
offer the best explanation for the complex variability in L2 Pup. 



3.4 Pulsation mode 

Jura et al. (2002) cited the difference between 12 /im light curve 
and those at H and K for L2 Pup as evidence that the pulsations 
are non-radial. Smith (2003) also noted the intriguing difference 
between the shorter and longer wavelength infrared observations 
of the DIRBE instrument, which was exceptional in a sample of 
207 infrared sources. In particular, the maxima at 1 .25 fim preceded 
those at 4.9 /im by 10-20 days, while at 4.9 fim, there was a sec- 
ondary peak between the two 1.25 ^m maxima. Also, the 2.2 ^m 
and 3.5 /im light curves resembled the 1.25 /im curve, while the 
12 fim curve was similar to that at 4.9 fim. The Jura et al. hypoth- 
esis was based on the assumption that, while the 12 /im flux tracks 
the light emitted by the circumstellar dust, and hence measures the 
luminosity of the entire star, the near-infrared fluxes measure only 
emission from the hemisphere of the star that faces the Earth. If 
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there are differences then the photosphere must vary in time non- 
spherically, implying non-radial pulsation, which would also natu- 
rally explain the marked time variation of the position angle of the 
net polarization of the star (Magalhaes et al. 1986). 

However, as wc have shown in Paper I, the extinction- 
corrected K band magnitude of L2 Pup is very close to the value 
expected from the Mira P-L relations. The good agreement tends 
to argue against the hypothetic nonradial pulsations (unless one al- 
lows that Miras also pulsate non-radially). The full velocity ainpli- 
tude of 12 km s^^ (Lebzelter et al. 2005) is also hardly compat- 
ible with non-radial oscillations, which are not expected to result 
in large radial velocity amplitudes (Wood et al. 2004). Looking at 
DIRBE light curves of Mira stars, it is apparent that not only are 
there similar phase lags in Mira stars (Smith et al. 2002), but also 
that small differences do exist between the near- and mid-infrared 
data. As discussed by Smith et al. (2002), there are no theoretical 
models that predict these phase lags, because of the complexity of 
the problem. In addition to the pulsation-induced luminosity varia- 
tions, the 12 jim light curve also reflects the dynamics of dust for- 
mation process which is governed by a time scale of its own (Win- 
ters et al. 2000 and references therein). Furthermore, the DIRBE 
observations were taken a few pulsation cycles before the begin- 
ning of the major dimming event in 1994. It seems to be likely that 
at the time of the DIRBE observations, there had already been some 
activity in dust formation, which can explain the extra features in 
the mid-infrared data without invoking non-radial pulsations. For 
example, Soker (2000) examined the effects of magnetic cool spots, 
above which large quantities of dust are expected to form. Alterna- 
tively, interactions with a close binary component (but outside the 
AGB envelope) may yield a strong density contrast in the equatorial 
plane, with similar outcomes. While currently there is no evidence 
for binarity in L2 Pup, Ireland et al. (2004), using optical inter- 
feromety, found evidence for highly clumpy structures in the cir- 
cumstellar dust shell as close as the dust condensation radius. Very 
recently. Wood et al. (2004) presented observations that suggested 
large scale star spot activity in some semiregulars with long sec- 
ondary periods. Although L2 Pup does not belong to that subclass 
of SRs, we think that the Soker (2000) model offers a more realistic 
mechanism to produce inhomogeneous dust and therefore different 
mid-infrared light curves in L2 Pup than non-radial pulsations. 



4 CONCLUSIONS 

L2 Pup is a very interesting semiregular star. Its light curve is dom- 
inated by two completely different mechanisms: gradual dimmings 
caused by circumstellar dust and pulsations. These two are com- 
pletely independent, as we did not find any evidence for coupling 
between them. The present paper discussed properties of the pulsa- 
tions, outlining the physical implications. 

The frequency spectrum of the long-term light curve shows 
a peak of power that is resolved into multiple peaks under a nar- 
row envelope. This is consistent with stochastic excitation, as seen 
in solar oscillations. The mode life-time can be derived from the 
width of the envelope: for L2 Pup this is found to be about 5 yr. 
The oscillation amplitude agrees roughly with the predictions of 
simple scaling laws, while the phase behaviour is markedly differ- 
ent from a Mira star with the same pulsation period. The seemingly 
random amplitude and phase changes argue against excitation by 
the K-mechanism alone, although it might make some contribution. 

A test for low-dimensional chaos ruled out the possibility of 
non-linear interactions of a few pulsation modes as the reason for 



seemingly irregular light variations. We also argued against the hy- 
pothesis of non-radial pulsations in L2 Pup: the difference between 
the near- and mid-infrared data may only be related to clumpy dust 
production, possibly driven by spot activity. 

Bedding (2003) has recently given examples of other semireg- 
ulars whose power spectra display similar evidence for solar-like 
oscillations. The inferred mode lifetimes range from a few years, 
as found in L2 Pup, down to less than a year. It seems plausible 
that stochastic excitation, presumably from convection, plays an 
important - perhaps dominant - role in the behaviour of semiregu- 
lar variables. 
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